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All beams of electromagnetic radiation are made of photons. Therefore, it is important to find 
] a precise relationship between the classical properties of the beam and the quantum characteristics 

(<~*) , of the photons that make a particular beam. It is shown that this relationship is best expressed 

' in terms of the Riemann-Silberstein vector — a complex combination of the electric and magnetic 

' field vectors — that plays the role of the photon wave function. The Whittaker representation of 

this vector in terms of a single complex function satisfying the wave equation greatly simplifies the 
analysis. Bessel beams, exact Laguerre-Gauss beams, and other related beams of electromagnetic 
radiation can be described in a unified fashion. The appropriate photon quantum numbers for these 
beams are identified. Special emphasis is put on the angular momentum of a single photon and its 
connection with the angular momentum of the beam. 
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"The notion of photon as a quantum of electromagnetic energy, like the notion of atom as an elementary unit of 
matter, permits a number of interpretations. " With these words Bruce Shore begins the section entitled "What is a 
photon" in his comprehensive monograph and then he goes on to present a few possible answers to this question. 
We hope that our contribution to this special issue adds a little to the neverending discussion on the nature of photons. 

Electromagnetic radiation, especially in the optical range, is produced and used most often in the form of beams. 
There is a variety of mathematical models to describe such beams: the Bessel beams, the Hermite-Gauss beams, the 
IL^ Laguerre-Gauss beams, and also the focus wave modes. The mathematical representations of these beams are the 
exact, or approximate, solutions of the classical Maxwell equations. In recent years, there were many experiments that 
exhibited the influence of the orbital momentum on the beam properties (phase dislocations, optical vortices) and on 
?H ■ the interaction of beams with matter (particle trapping, optical tweezers and spanners). As a rule, orbital angular 
momentum leads to vortices and electromagnetic beams with vortices, as we shown recently |^|^, are interesting 
because they can guide charged particles. The fundamental theoretical and experimental papers on the optical angular 
momentum were recently reprinted in a collection In all these papers, the theoretical description of beams carrying 
angular momentum is given in terms of the classical solutions of the Maxwell equations and sometimes separately in 
terms of photons. However, no unifying principle is given that would connect these two points of view in a precise 
manner. The purpose of this work is to fill this gap. We shall show that the notion of the photon wave function 
appears in both of these descriptions and that it provides a very convenient concept to unify the t wo p oints of view. 
A very useful mathematical tool in this analysis is the Riemann-Silberstein vector |a, Ifl la. S llU UM ■ Applications 
of the RS vector to many physical problems were recently reviewed in a very thorough paper by Keller jl2| |. The 
Whittaker representation [l^ of this vector greatly simplifies the calculations since then the vector field is described 
by a single function. We shall not consider here the beams described by approximate solutions of Maxwell equations 
obtained in the paraxial approximation p^ . The analysis of such solutions in terms of photons would be rather 
awkward — the notion of an approximate photon does not make sense. 



'Electronic address; |birula@cft.edu.pTl 



2 



There is some overlap between our final conclusions and those obtained recently by Jaregui and Hacyan jlSj even 
though our methods are completely different. They rely on the standard description of the quantized electromagnetic 
field based on the vector potential. We follow the methods developed in our earlier works on quantum electrodynamics 
and photon wave functions [TH where the Riemann-Silberstein vector plays the central role. 
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II. SUCCINCT DESCRIPTION OF THE ELECTROMAGNETIC FIELD 

A natural tool in the analysis of the solutions of the Maxwell equations, in both classical and quantum theories, is 
the Riemann-Silberstein (RS) vector F 

^{E + icB). (1) 

The physical significance of the RS vector has been recognized by Silberstein ^ who observed that important char- 
acteristics of the electromagnetic field (energy density, Poynting vector, Maxwell stress tensor) are bilinear products 
of the components of this vector. The total energy 

E = j (frF*-F, (2) 

the total momentum (i.e. the integral of the Poynting vector divided by c) 

£rF*xF, (3) 

J 

and the total angular momentum 

M = — Jd^r {rx {F* x F j) (4) 

look very much like the quantum-mechanical expectation values and we shall show later that this fact has a deeper 
meaning. 

The convenience of using the RS vector has been recognized by Bateman who was the first to analyze with its 
help various solutions of the Maxwell equations. Kramers 16\ used the RS vector to formulate the canonical theory of 
the electromagnetic field and Power [13 stressed the usefulness of this vector in the description of circularly polarized 
waves. The complex RS vector carries exactly the same information as two real field vectors but its use significantly 
simplifies the mathematical analysis. In particular, the two pairs of real Maxwell equations written in terms of F 
reduce to just one pair of complex Maxwell equations 

dtF{r,t) = -icV X F{r,t), V-F(r,t) = 0. (5) 

The RS vector can be expressed in the following form 

F(r, t) - V X (ldtZ{r, t)+Wx Z{r, t)^ , (6) 

where Z{r,t) is a complex vector field (a unified form of the two Hertz vector potentials) satisfying the d'Alembert 
equation 

(la,2-A)Z(r,t)=0. (7) 

In the description of beams, it is convenient to choose the vector Z in the direction of propagation Z{r, t) = 
(0, 0, 1)x(t, t)). In this way, we obtain the following representation of the RS vector in terms of one complex function 

= {d^d,+ '-dydt)x{r,t), (8a) 

Fy = {dyd,-'-d,dt)x{r,t), (8b) 
= ^{dl + dl)x{r.t). (8c) 
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We shall refer to these formulas as the Whittaker representation. A century ago Whittaker discovered that an 
electromagnetic field obeying the Maxwell equation can be described by two real functions. By separating Eqs. © 
into the real and imaginary part, we recover the original Whittaker's formulas. The Whittaker representation (jHl of 
the RS vector is the simplest but it is not unique. For example, by choosing the complex Hertz vector Z in the form 
Z{r,t) = {l,i,0)x{r,t)) we obtain a different expression for F that will be used in Section Hxl 

The representation of the solutions of the Maxwell equations in terms of a single complex function (or two real 
functions for that matter) satisfying the d'Alembert equation greatly simplifies the analysis. This has been noticed 
already by Stratton 8] in his derivation of the formulas for the Bessel beams. Since there are no auxiliary conditions 
(like the vanishing of the divergence) imposed on x, this single function describes true degrees of freedom of the 
electromagnetic field. Sometimes, following Whittaker, we shall call x a "scalar" solution of the d'Alembert equation. 
However, strictly speaking this terminology is not justified because x has fairly complicated transformation properties, 
quite different from those of a scalar field. 

Every solution of the d'Alembert equation can be decomposed into plane waves. Therefore, we can write 

X{r,t) = jdkN{k) (/+(fc)e-'"^'=*+''="'+/-(fc)e*'^'=*-'''^"') , (9) 
where ujk = kc, N{k) is a normalization factor, and 

""-^^ (10) 

With the use of Eqs. (jSJ, we obtain the following decomposition of the RS vector into plane waves 

F{r,t) ^ Jdke{k) (/+(fc)e~*'^''*+*'' + /-(fc)e"^'=*-*'=-'') , (11) 
where the complex polarization vector 

kxkz + ikky \ 

-kykz ~ ikkx (12) 

^x+^y J 

is a normalized solution of the following set of algebraic equations 

n X e{k) — —ie{k), n — k/k. (13) 
We choose the normalization factor N{k) as 

1 



N{k)^^=^, k^ = ^kl + kl (14) 

to make e(fc) normalized to one (e(fe) • e*{k) = 1). Eqs. (|13|) together with the normalization condition define the 
polarization vector up to an arbitrary fc-dependent phase. This reflects the gauge freedom of the complex Hertz 
potential Z{r, t). 

The two arbitrary complex amplitudes f^{k) describe the true degrees of freedom of the electromagnetic field. 
They will be labeled with the index a taking on the values ±1. The sign of a determines whether the pair of vectors 
{E, cB) obtained from the plane wave solutions e(A;)e~*'^'^'^*~*'^'''-' of Maxwell equation (|3J) rotates clockwise (a = 1) 
or anticlockwise (cr = -1) along the direction of propagation. Thus, the amplitudes f^{k) and f~{k) describe two 
circular polarization components of the wave. 



III. BESSEL BEAMS 



In the previous section we have given the solutions of the Maxwell equations in the plane wave basis. Since we are 
interested in the beams with angular momentum, we shall expand the plane waves appearing in (jSJ according to the 
formula (Eq. 8.511.4 of Ref.|19j) 

OO 

^^k■r ^ ^ik^z i^e''^^'^-f^Jr,^{k^p), (15) 

m— — OO 
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where we introduced a convention that <f) is the polar angle in the xy plane and is the polar angle in the k^ky 
plane. With the use of this formula, we obtain an expansion of x('", i) into the solutions of the d'Alembert equation 
in cylindrical coordinates 

X{r,t)^ ^ / _£ / ^ (Xfc^mfc. (P^ '/'^ t)fk^,nk, + Xfc^mfc, (P. '/'^ t)fk^,nk,)^ (16) 

where 

Xlk^rniP, ^, ^, t) = ^=^e — J™(fc^p), (17) 



and 



2ir 



= d^e— '^r(fe). (18) 



The RS vector, calculated from Eqs. ((HI with x given by (|17|l . has the form 

{^a) 



ifTfc_(a)e-*J™+i(e)+iafc+(a)e-'^^J™_i(?) 
F^^„feJp,</.,z,i) = Si^e-'^("'=*-'=^^-"« ( fc_(a)e-^J™+i(e)-fc+(a)e— ^J„_i(0 |, (19) 



where = k±p and 



k±{a) = (20) 



In the derivation of p9|l we have used the following recurrence relations for the Bessel functions 

2mJ„(e) = Un-iiO + ^Jm+iiO, (21a) 

2d^JmiO = Jrn-liO~Jm+liO- (21b) 

The Bessel beams p9|l form a complete set — all other solutions can be obtained as linear combinations of Bessel 
beams. The Bessel beams are characterized by the four parameters kz,k±,m, and cr. The parameters kz and fcj^ 
{—OO < kz < OO, < k± < oo) are the z and the p components of the wave vector, respectively. The parameter 
m takes on all integer values. The components of F^^^-^j. in the cylindrical coordinates are simpler and allow for a 
complete factoring of the angular dependence 




/. / iakzd^+ikm/^ 

V ) I -(jfcS^-fc^m/e I J™(0- (22) 



V2k 



Note, that the phase factor in front is invariant under the following simultaneous changes of the coordinates z and 0: 
z —^ z — ma/kz, (f) ^ (f) + a. That means that the Bessel beam has a screw symmetry. 

Bessel beams take on a very simple form in the limit, when k^ —> 0. In order to avoid the trivial (zero) result, we 
shall divide (for m > 0) the expression (|17l) by fc™~^ before taking the limit 

Yiu,{k\--\l^„,{pA,z,t) = ^ > e"-(l^'l-*-^--)(x + zay)"'. (23) 
fc-L^o V2fc2™m! 

This expression represents on the one hand an exact solution of the wave equation but on the other hand, when 
substituted into the Eqs. ||HJ), it also gives an approximation of the Bessel beam in the vicinity of the z axis. The 
electromagnetic field derived from H23|l has some unique properties j20]. In particular, when to = 2, it allows for exact 
solutions of the equations of motion for charged particles 



•ope 




IV. WAVE PACKETS OF BESSEL BEAMS 



Bessel beams are almost ideal representations of pencil-shaped electromagnetic radiation moving along the z axis 
without spreading in the transverse direction. In this respect, they are a much better representation of realistic 
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electromagnetic beams than plane waves. However, there is a price to be paid for the diffraction-free propagation: the 
fall-off of the field is so slow (like that the energy of the electromagnetic field per unit z is infinite. A standard 

method of avoiding the infinite energy problem for monochromatic beams consists of abandoning the requirement 
that the functions describing the beams are exact solutions of the Maxwell equations. Along these lines, various 
forms of Gaussian beams, i.e. beams with a Gaussian fall-off in the transverse direction were proposed (see We 
find the departure from Maxwell equations (paraxial approximation) to be too high a price to pay for saving the 
monochromaticity — strictly monochromatic beams do not exist in reality, anyway. We shall show that a fast fall-off 
in the transverse direction can be obtained for exact solutions of Maxwell equations constructed as wave packets of 
Bessel beams. However, these beams will be, of course, not strictly monochromatic. 

Instead of working with the RS vectors, we shall form the wave packets using the scalar functions x- Since the 
transformation (jS)) from % to is linear, the wave packets for the complete RS vector can be easily produced by taking 
the appropriate derivatives of x- Taking this into account, we shall often use the name "beams" for scalar solutions of 
the d'Alembert equation with the understanding that the solutions of Maxwell equations describing electromagnetic 
beams are to be obtained from x according to Eqs. 

We shall consider the most general wave packet of Bessel beams with given values of m and a. A wave packet 
solution of the d'Alembert equation of this type can be written as the following superposition of the functions l|17|l 



Xm.{pA,z,t)^e'^'^^ dk^k^ dfc,5(fc„fci)e-'"(V'=Hfci*-fc^-)j„(fc^p), (24) 



where g{kz,k±) is some weight function. For beam-like solutions, the longitudinal component k^ should be much 
larger than the transverse component k±. Note, that this condition is also essential for the paraxial approximation. 
In the limit, when the weight function shrinks to a product of delta functions 

g{k.,k^)^^j^J{K.-k'Mk±-k'^). (25) 

we obtain back the function (|17|l . There are not too many weight functions that allow for an explicit evaluation of the 
integral (|24|l but there are several examples when this can be done. Not all of them will lead to well formed beams. 
The most important case is treated in the next section. 



V. EXACT LAGUERRE-GAUSS (LG) BEAMS 

Exact Laguerre-Gaussian beams are most easily obtained by changing the variables in the formula (|24|l from k^ 
and fc_L to k± = {lOhjc ± kz)l2 = (fc ± kz)l'2.. We shall also introduce the new variables t± = t:k z/c and rewrite the 
exponent in Eq. H24|l in the form 

uJkt — kzZ — c{k+t^ + k^t+). (26) 
In these variables the general wave packet with a given values of m and a has the form 

X;;(p,0,z,t) = e^-"M / dfc+dfc_g(fc+,fc_)e-^"='=+*-e-^'^^^--*+J™(2yE^p). (27) 
Jo Jo 

The dependence on t-(i+) signifies the propagation of the signal in the positive (negative) direction of the z axis. In 
a general superposition of Bessel beams ()27|l . both these directions of propagation are present. However, it is clear 
that in a true beam-like situation one of the propagation directions must play a predominant role. We shall choose 
the convention that the beam propagates in the positive z direction. In this case, for a true beam the dependence 
on must be weak and only cause a modulation of a propagating wave. In the limiting case, when there is no 
dependence on tj^^ the wave propagates along the z axis without changing its shape. In order to construct a beam- like 
wave packet, with fc+ concentrated near some large value O/c, we shall allow only for a narrow range of values of fc_ 
as compared to Vt/c. 

The simplest functions g{k+, fc_) that generate exact beam-like solutions of the Maxwell equations with a Gaussian 
fall-off in the transverse direction have the form 

g{k+, fc_) = 6{k+ - r!/c)fc"+™/'e-''"'=-/'; (28) 

where / , as will be shown below, determines the width of the LG beam. With this choice of g (fc_|_ , fc_ ) , both integrations 
in Eq. (|27|l can be performed. Owing to a delta function, the integration over k-^ is immediate and it results in the 
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replacement /c+ fl/c. Next, we shall integrate over fc_ using (Eq. 6.643.4 of Ref.[ll) 

I d..-^^l^e-^M2PV-x) - K (^), (29) 



where L™ is the associated Laguerre polynomial. This gives 

where a(t+) = + icrc^t+Zfl and A = n!(c/rj)"+'"/2+i. Finally, our wave packet of Bessel beams representing an 
exact LG beam, expressed in Cartesian coordinates, takes on the form 

X^nnn.{p, c^, z, t) = A ^(^^)„+,„+i exp (^"^ j j • (31) 

The width of the Gaussian is equal to I at the waist and grows with as we move away from — 0. In the 
optical regime, when = lO^^sec"^ and I = 0.001m, the modulation of the beam through its dependence on t+ has 
characteristic time l^n/c^ « lO^^sec^^. This time is seven orders of magnitude longer than the wave period. 
In the simplest case, when n = 0, m ~ 0, the exact LG beam reduces to a pure Gaussian 

X^nooi^,y,z,t) = ,(,^) exp (^-^j . (32) 

Apart from some changes in notation, this function for a — 1 coincides with the solution investigated recently by 
Saari, Menert, and Valtna in connection with photon localization. Their solution, as it turns out, is a member of 
a complete set of functions, labeled by n and m, all having an exponential fall-off in the transverse direction. There 
is a vast literature on this subject and various forms of the solutions of the wave equation equivalent to our formulas 
can be found in What makes our approach different is a systematic use of the Bessel beams as a basis from 

which other solutions can be obtained. 

Obviously, for each function Xcinm obtain from the formulas JHJ a solution of the full Maxwell equations. The 
exact LG beams fall off sufficiently fast in the transverse direction to guarantee that all physical quantities (energy, 
momentum, angular momentum) per unit interval in the z direction are finite. 



VI. SPECTRAL DECOMPOSITION OF EXACT LAGUERRE-GAUSS BEAMS 

The exact Laguerre-Gauss beams given by the formula H31|l are similar to the so called elegant LG beams |23ll2^ . 
However, in contrast to the elegant LG beams, the exact LG beams are not monochromatic but they solve the 
Maxwell equations exactly and not only in the paraxial approximation. Even though our exact LG beams are 
not monochromatic, they can be made nearly monochromatic by a proper choice of the parameters f2 and I — their 
spectrum will be sharply peaked. In order to see this, we shall decompose the function H31|l into its Fourier components 
in the time variable 

1 /"^ 

Xlanmi^^y^z) = — j dte"^\'^„^{x,y,z,t). (33) 

This task is made simple by returning to the original formula (|27|) and observing that the integration over time leads 
to a delta function (5(c(fc+ + fc_) — auj). The sign of lo must coincide, therefore, with the sign of cr, i.e. a = sign(ti;). 
The absolute value of uj is never less than il because k+ = fl/c and A:_ > 0. The final expression for the Fourier 
transform, obtained after using the two delta functions, can be written as a product of the spectral weight function 
w{lo) and the time-independent part of the Bessel beam 

X^nn^iP, <l>, z) = w;(c^)e*'=^^/=e""^ J,„(fcxp), (34) 
where the spectral weight function is 
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FIG. 1: The normalized spectral weight w{lj) plotted as a function of uj for the following values of the parameters: a — 1, 
= lO^^s"^ and I = 0.001m. The sharpest peak represents the case n = 0,m = and the remaining two correspond to 
n = l,m = 1 and n = 2,m — 2. As we increase these numbers, the peaks flatten more and more but they are always localized 
near Q. 



C is a constant, and the parameters kz, kj_ and fc_ arc now the following functions of Q and to 



fc^ = (2f7- |a;|)/c, ki_ 2y/{\Lu\ - n)n/c, k- ^ {\uj\ - n)/c. (36) 
In order to make the comparison of different cases easier, we will normalize all weight functions to 1 

J dLow{u}) = 1, (37) 

which changes the value of the constant C to 

^ " r(n + m/2 + 1) • ^^^^ 

The spectral weight function of an exact LG beam has a cutoff at = f7 and is concentrated near this frequency. 
The prefactor /j""*"™/^ shifts the maximum from il to + c(n + to/2)/2/. In Fig. ^ we show the normalized weight 
functions for three sets of numbers n and m. 



VII. QUANTIZED ELECTROMAGNETIC FIELD 



In quantum electrodynamics the RS vector and its complex conjugate are replaced by the field operators built from 
the operators of the electric and magnetic field vectors 



F{r,t) 



■icB{r,t) 



(39) 



The canonical quantization of the electromagnetic field that employs the mode expansion into a set of harmonic 
oscillators |2^ |2g yields the following equal-time commutation relation between these field operators 



F,{r,t),FUr',t)\ = -hce,jkdk5^^\r - r'), \F,{r,t),Fj{r',t) 



-0, 



Fl{r,t),F]{r',t) 



= 0. 



(40) 



The expansion of the RS operator into plane waves is obtained from Eq. Hll|) by replacing the classical amplitudes 
/^(fc) by the creation and annihilation operators f^{k) —> \/hd{k) and ,f^{k) \/hb^{k) 



F{r,t) = Vh jdke{k) (a{k) 



^ — iLot-\-ik-r 



ot—ih-r 



(41) 



The amplitude / (fc) is being replaced by the creation operator because it is multiplied by the opposite frequency 
factor exp{iu}t). Commutation relations between the creation and annihilation operators 



[a(fc),at(fc')] ^{'2nfuj5^^^\k-k'), \b{k),h\k')\ ^ {2nfuj 5'^'^\k - k'), 



(42) 
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with all remaining commutators vanishing, reproduce correctly Eqs. H40() since 

ei{k)e*{k) - e*{k)ej{k) = ieijknk- (43) 



The creation and annihilation operators in this work differ by a factor {2tt)^uj from the traditional ones that satisfy 
the commutation relation (|42|) with just a straight delta function on the right hand side. We have departed from the 
standard normalization of these operators for two reasons. First, we wanted to keep a close connection between the 
classical amplitudes and the annihilation/creation operators. Second, as shown in [9j, our normalization makes the 
photon density operator in momentum space a^(fc)a(fc) a relativistic scalar. 

The commutation relations H42(l enforce the replacement of the classical amplitude f^{k), that multiplies the 
negative frequency solution, by the creation operator w{k). Replacing it by the annihilation operator would yield the 
wrong sign in the commutation relations H42() . The operators a^(A;) and fo^(fc) create different photons (of opposite 
circular polarization). The fact that the RS field operator contains only a{k) and b^{k) (and not b{k) and a'l'(fc)) 
leads to some unique properties of this field operator which are different from those of the electric and magnetic field 
operators. To see this, let us consider a general state of a single photon, created from the vacuum by an arbitrary 
linear combination of creation operators and 

«l = /— h^{k)aHk) + 4,-{k)b^k)) . (44) 
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We have to impose the normalization condition 

(|^+(fc)|2 + |^-(fc)|2) ^1 (45) 

to guarantee proper commutation relations between and . With the use of this creation operator we may define 
a coherent state of the electromagnetic field 

l^)c - cxp (-^(^>) exp |0), (46) 

where (N) is the average number of photons in the coherent state. The expectation value of the field operator in a 
coherent state is given by the formula 

(F{r,t))^ = ^/h{^) Jd\e{k){^+{k)e~"^'+'''-''+{i;-{k)ye''^'~'''-'') . (47) 

Note the complex conjugation of the wave function in the second term. This is caused by the fact that this term 
has opposite sign of the frequency (negative energy in the quantum- mechanical setting). Every mode function in the 
expansion of the RS vector with negative energy must be interpreted as the complex conjugate of the photon wave 
function. 

Since coherent states are the right eigenstates of the annihilation operators and the left eigenstates of the creation 
operators 

a(fc)|V)c = ^+(fe)|^)c, c(V'l&^(fc)=c(V^|V'"*(fe), (48) 

the expectation value of any product of the RS field operators is equal to the product of the corresponding classical 
RS vectors. One can show with the use of that the expectation value of any departure from the average value 
vanishes 

((F(r,<)-(F(r,i))e)")c = 0. (49) 

Note, that this property holds just for the RS operator (and also for its hermitian conjugate). It is so, because the 
RS operator contains the annihilation operators of right-handed photons and the creation operators of left-handed 
photons. The electric and magnetic field vectors exhibit always (even for coherent states) large fluctuations around 
the average value. 

The formula (|47|l gives also an inverse connection between a classical RS vector and a quantum coherent state. 
From a solution of classical Maxwell equations, we may find the coefficients in an expansion of the RS vector into plane 
waves. Their norm determines the average photon number (TV) while the normalized coefficients ip±ik) determine 
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the creation operator (|44() . Together, they enable us to define the coherent state (|46ll . Obviously, the average value 
of the RS field operator in this coherent state reproduces the original solution of the Maxwell equations. 

Since in this paper we are interested in the properties of beams carrying angular momentum, we shall now use the 
expansion of the operator F{r,t) into the Bessel beams. Following the same steps as in the classical case, we obtain 

F{r,t) = VhY^ / / ^-i(F,+ ,„,^(p,,^,z,t)a(/ci,m,fc,) + F,7^„,Jp,0,z,t)6t(fc^,m,fe,)), (50) 



m— — oo 



where the new creation and annihilation operators are defined as 

a){k^,m,k,)^^ I dipe"'"^a''{k), a{k±,m,k,) = ^ [ dip e''""^ a{k) , (51a) 
27r Jo 27r Jo 

b\k^,m,k,) = — d^e'^n^k), h{k^,m,k,,) = — d^e-'^n{k). (51b) 

Thus, F^^^f, play the role of mode functions. According to the standard interpretation of field operators in relativistic 
quantum field theory, -F)!^^^,, is the wave function of the photon annihilated by d{k±,m,kz) and -F/"^^^.^ is 

the complex conjugate of the wave fimction '4'kj_mk ^^"^ photon created by &^(fc^, m, kz)- 

.m ( ^fc_e** J™+i(0 + ik+e-'"^ J,n-i{^) \ 

( -ik+e^'^Jm+iiO - *fc-e-'*J,„_i(e) \ 



k±Jm{£,) 



The notion of the photon wave function has been discussed in detail in our review paper s-nd also, more recently, in 
the report by Keller [T^ . In the next section we shall analyze the quantum mechanical properties of the single-photon 
states described by the wave functions V'fc'^mfe^- 

VIII. PROPERTIES OF PHOTON WAVE FUNCTIONS 

We have constructed in the previous section photon wave functions using the apparatus of quantum field theory. 
However, there is a more pedestrian way that may be useful to better explain this construction. Let us assume that 
we want to repeat the success of Schrodinger and we want to construct a wave mechanics of photons. We shall need 
photon wave functions xp satisfying an analog of the Schrodinger equation and a set of operators representing basic 
physical quantities. Since photons have spin one, we need a vector wave function. The general form of the Schrodinger 
equation is 

ihdtxp{r,t) = H'ip{r,t) (54) 

and all we need is a photon Hamiltonian. The operators representing physical quantities associated with the gen- 
erators of basic transformations are easily constructed. The momentum operator (associated with an infinitesimal 
displacement) and the angular momentum (associated with an infinitesimal rotation) have the form 

p=-V, M^-rxV + hs, (55) 
i i 

where s stands for a vector built from the spin-one matrices, 

i\ / -i \ 

(56) 

We shall also need the helicity operator A — the sign of the projection of the angular momentum on the momentum. 

A= , \ ={s-p). (57) 

^pI+pI +pI 
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It has been noticed seventy five years ago (cf. Ref.^Jl for a historical review) and rediscovered over and over again, 
that the complex form of Maxwell equations © with the help of the relation 

Vx=-i{s-V) (58) 

can be cast into a Schrodinger-like form 

ihdtrpir, t) = c{s-p)ip{r, t), (59) 

with the rescaled helicity operator playing the role of a Hamiltonian H — c{s-p). Thus, there is a close relation 
between the helicity and the Hamiltonian. However, one must exercise care because the Maxwell equations have 
solutions with positive and negative frequencies, while photons always have positive energy. Therefore, solutions of 
the complex Maxwell equations with negative frequencies must be interpreted as complex conjugate wave functions. 

We shall now analyze the wave functions M't^^mk associated with the Bessel beams using the tools provided by wave 
mechanics. To this end we choose the following set of four commuting operators: the component of the momentum 
in the z direction p^, the square of the transverse momentum the component of the angular momentum in 

the z direction M^, and the helicity A. It is only a matter of straightforward calculations to check that both wave 
functions V'fej^mfe^ are eigenfunctions of the first three operators belonging to the eigenvalues hk^, h^k^, and hm, 
respectively. To verify that they are also eigenfunctions of the helicity operator, no additional calculations are needed. 
Both functions ipk±mk satisfy the complex form of Maxwell equations, but they have opposite signs of the frequency 
ujk- Therefore, we obtain V x i/)^ = ±kip^ , and since there is no imaginary unit in these equations, the same relation 
holds for complex conjugate functions. With the use of l|58|) . we obtain 



^'^k I mk, = T^/^k^mk^ — V'fc , mk^ = ^'^k , mk, ■ (^0) 



Thus, according to the discussion in the previous section, a classical Bessel beam (|19|l is made of photons with the 
quantum numbers hk^, h^k\, hm, and cr. 

There is one problem with the photon wave function that for many physicists invalidates this notion altogether. 
Namely, the standard Born interpretation of the photon wave function xp{r^t) fails because there is no well defined 
position operator. This problem, however, is common to all particles described by relativistic quantum mechanics 
and it simply follows from the fact that a positive energy solution ip^{r, t) of any relativistic wave equation loses this 
property when it is multiplied by r. In particular, we encounter this problem for electrons described by the Dirac 
equation. In this case, however, the probabilistic interpretation is based on the charge distribution. Since photons do 
not carry any charge, we have to base the probabilistic interpretation on the energy distribution. 

The RS vectors can serve as photon wave functions to calculate an expectation values of a quantum-mechanical 
operator O provided this expectation value is defined in the following way 0, 

(6) = iV-i / (frF*{r,t) ■ —^d^=F{r,t). (61) 

1^1 J\H\ 



where the absolute value takes care of the fact that the Hamiltonian has both positive and negative eigenvalues and 
the norm N is defined as 

N ^ j (f'rF*{r,t) ■ ^F{r,t). (62) 

This prescription can be intuitively understood as follows: The energy density F*-F must be "divided" by the energy 
to produce a "probability density". 
The division by the energy can be implemented in Fourier space as a division by hujk- The resulting expression, 
when transformed to coordinate space has the form 

^-tA^ fd'r [d\'F*{r,t)-r-^F{r',t). (63) 



27r . . 

Using the formula H61fl one can check that the correspondence principle is satisfied for the energy, momentum, and 
angular momentum, namely 

(H) = 1 jd\F*-F, (64a) 
(p) = {-ihV) = ^ J(frF* X F, (64b) 
(M) = {-ihr xV + hs):^^ I d^r {r X {F* X F)) . (64c) 

Ci V J 
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The formula for the expectation value of the helicity, like the one for the norm, has a nonlocal character 

(A) = (^)-^/^^^^VF*(r,.).^VxFKt). (65) 

This double integral has been obtained by Deser and Teiteltboim 28] as a generator of dual rotations in classical 
electrodynamics. 

Let us suppose now that F is an arbitrary eigenfunction of Mz ~ —ih{xdy — ydx) + fiSz belonging to an eigenvalue 
hm. In this case, from Eq. H64c() we obtain 

{M^) = — Jd^r {rx {F* x F))^ = Nhm. (66) 

The number N, playing the role of the norm in quantum mechanics of the photon, acquires a physical interpretation 
in quantum electrodynamics Isol Isif . Namely, it is the average number of photons N = (N) in a coherent state 
associated with the RS vector F. Hence, the value of the angular momentum obtained from the classical formula is 
equal to the quantum eigenvalue hm (i.e. the value per one photon) multiplied by the number of photons. 



IX. BEAMS IN MOMENTUM REPRESENTATION 



Momentum representation, well known from nonrelativistic quantum mechanics, is very useful in the description of 
photons since it is free of the difficulties with the probabilistic interpretation encountered in the coordinate represen- 
tation. Actually, we have already introduced the momentum representation by writing down the formula pi|l for the 
RS vector. The amplitudes if^'^ik) — f^{k) and V' (fc) = (/ (fc))* may be viewed as the photon wave functions in 
momentum representation. Complex conjugation in the second formula, as we have explained in Section [VIII is due to 
the fact that the amplitude f~{k) is multiplied by the plane wave factor with opposite frequency. The amplitudes ^/^^ 
carry full information about the RS vector. Every operator O acting on the left hand side of Eq. (|ll|l can be translated 
into an operator in the momentum representation acting on the amplitudes ?/'^(fc). In particular, an infinitesimal dis- 
placement V acting on F{r, t) is reproduced by a multiplication of ?/;^(fe) by ik. Therefore, the momentum operator 
p acts on ^^{k) simply as a multiplication by hk, as is always the case in momentum representation. The action of 
an infinitesimal rotation r x V -I- is on F{r, t) translated into an action on ip^lk) is not unique. Its form depends on 
the choice of the polarization vector. For our choice of e(fc), that resulted from the Whittaker representation ||SJ), we 
obtain the following expressions for the components of the angular momentum operator 

M,V^(fe) = h(^-i{kxdkh±^^^^i;^{k), (67a) 

/ kk \ 

M,^±(fc) = h -i{k X dk)y ± j^-^ ^±(fc), (67b) 



1,2 4_ 1.2 



Mz^^{k) = -ih{kx dk)zi}^{k). (67c) 

The Whittaker representation is very convenient because in this representation Mz reduces just to the z component of 
the ordinary orbital angular momentum so that the eigenstates of Mz depend on the polar angle f through the phase 
factor exp(zTO(^). This may look counterintuitive because one would expect also a contribution from the photon spin. 
The absence of such a contribution is a consequence of the special choice of the phase of e(fe) and hence of the phase 
of /^(fc). Another natural choice of the complex Hertz potential mentioned in Sectional leads to the formulas 



M,^*(fc) = h^-i{kxdk)x±^^ji^^{k), (68a) 

My^^{k) = h(^-i{kxdk)y±^^^i^^{k), (68b) 

Mz^^ik) = h{-i{kx dk)z±l)^^{k). (68c) 

Here, the form of Mz agrees with our expectations. There is a contribution from the photon spin equal to ±h, but it 
was manufactured by choosing a particular phase of e(fc). 

In summary, only the total angular momentum of the photon has a physical meaning (as a generator of rotations) 
while the separation into the orbital and the spin part is not unique. The projection of the spin part on the direction 
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of momentum — helicity — (orbital angular momentum does not contribute here) is also well defined. For all choices 
of e(A;), the action of the helicity operator in momentum representation, namely A = k-M/hk, results in Atp^ = ±tp^. 

Now, we apply these results to the description of the Bessel beams and exact LG beams. In order to make a 
distinction between the labels and the arguments, we introduce a primed vector k' to denote the arguments (integration 
variables in the plane- wave representation) of the functions We retain the vector k for labeling the quantum 
numbers of a beam. For the Bessel beams we obtain the amplitudes f^{k') in the form 

ftmkSk'r^^'X) = (-*)"V2(2^)2fcV™'^'<5(fcl-fe^)5(fc;-fc.), (69a) 

fk^^mkSk'^^^'X) - {ftr..kSk'^^v'X)r- (69b) 

Indeed, upon the substitution of H()9al) into the formula (jllf) and after evaluating the integrals with respect to kz and 
fcj^, we can use the standard integral representation of the Bessel function to obtain 

^^^^{mv+k^pcos{v~<|>)) = 2^"" e'""*' J„,{k^ p) . (70) 







The formula for a = — 1 is obtained by complex conjugation. The amplitudes for the exact LG beam can be 
obtained by convoluting the expressions (|69|l with (|28|l and they have the form 

fLr,M^^',k'_) - (-^)'"^/2(2^)2fc'e™^'<5(fc;-^^/c)fc:"+'"/^e-'^"'='/^ (71a) 

Since the mode functions for the Bessel beams (|69|l and Laguerre-Gauss beams (I71|l describing opposite polarizations 
are related by complex conjugation, the corresponding wave functions are identical, as was to be expected. 



X. CONCLUSIONS 



We have built exact solutions of Maxwell equations carrying angular momentum as superpositions of Bessel beams. 
In particular, we analyzed the spectral properties of the exact Laguerre-Gauss beams. We have shown that we may 
associate photon wave functions in the form of the Riemann-Silberstein vector with all classical solutions of Maxwell 
equations. The Whittaker representation of the RS vector, that enables one to condense full solutions of the Maxwell 
equations into just one complex function, turned out to be particularly useful in this context. The quantum numbers 
of the photon states may serve as labels characterizing the classical solutions. The photon states with a prescribed 
angular momentum, such as the Bessel and Laguerre-Gauss beams, replace the commonly used plane waves and are 
best suited to describe light beams. The photon wave functions also appear in quantum electrodynamics as the mode 
functions in the expansion of the electromagnetic field operator into creation and annihilation operator. We have 
demonstrated that the RS field operator has a unique property: in coherent states the expectation value of the rt-th 
power of the RS operator is equal to the n-th power of the expectation value of this operator. This means that, in 
contrast to the electric and magnetic fields, the RS field operator exhibits no fluctuations around the coherent state 
average value. We identified the source of the difficulties in the splitting of the total angular momentum into the 
orbital and the spin parts. The nonuniqueness of this splitting is caused by an arbitrariness in the choice of the phase 
factor of the photon wave function in momentum representation. This, in turn, may be linked to the gauge freedom 
of the complex Hertz vector. 
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